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The part of the proton spin £ carried by u, d, s quarks is calculated in the frame- 
work of the QCD sum rules in the external fields. The operators up to dimension 
9 are accounted. An important contribution comes from the operator of dimension 
3, which in the limit of massless u, d, s quarks is equal to the derivative of QCD 
topological susceptibility x'(0)- The comparison with the experimental data on S 
gives x'(0) = (2.3 ± 0.6) X 10 -3 GeV 2 . The limits on £ and x'(0) are found from 
selfconsistency of the sum rule, £ = 0.33 ± 0.20, x'(0) = (2.4 ± 0.5) ■ 10~ 3 GeV. 
The values of gA = 1-37 ± 0.10 and = 0.65 ± 0.15 are also determined. 



Wc calculated EJ the value of £ - the part of nucleon spin carried by three 
flavours of light quarks £ = Au + Ad + As, where Au, Ad, As are the parts 
of nucleon spin carried by u, d, s quarks. On the basis of the operator product 
expansion (OPE) E is related to the proton matrix element of the flavour 
singlet axial current j° 5 

2ms^E = (p, s\j^ 5 \p, s), (1) 

where s^ is the proton spin 4-vector, m is the proton mass. The polarization 
operator 

n(p) = i J ^^(Oin^Cx), *7(0)}|0> (2) 

was considered, where T)(x) is the current with proton quantum numbers!. It 
is assumed that the term AL = j^A^ where is a constant singlet axial 
field, is added to QCD Lagrangian. In the weak axial field approximation H(p) 
has the form 

n(p) = n(%) + nW(p)v (3) 

n^ 1} (p) is calculated in QCD by OPE at p 2 < 0,\p 2 \ > R~ 2 , where R c is 
the confinement radius. On the other hand, using dispersion relation, njj (p) 
is represented by the contribution of the physical states, the lowest of which 
is the proton state. The contribution of excited states is approximated as a 
continuum and suppressed by the Borel transformation. The desired answer is 
obtained by equalling of these two representations. 



An essential ingredient of the method is the appearance of induced by the 
external field vacuum expectation values (v.e.v). The most important of them 
in the problem at hand is 

<0|j° 5 |0}^3/ 2 ^ (4) 

of dimension 3. The constant /q is related to QCD topological susceptibility. 
We can write 

(0\Ju 5 \0)a = Um q ^ i ( d 4 xe^(0\T{j%(x),j° (0)}\0)A V = 



= lim q ^ a P tlv (q)A v (5) 
The general structure of P^ u {q) is 

P„u{q) = -P L (q 2 )S^ + P T (q 2 )(S^q 2 + %q v ) (6) 

Because of anomaly there are no massless states in the spectrum of the sin- 
glet polarization operator P Mt/ even for massless quarks. Px^iq 2 ) also have 
no kinematical singularities at q 2 = . Therefore, the nonvanishing value 
Pfj,v(0) comes entirely from P^q 2 ). Multiplying P^(q) by q^q v , in the limit 
of massless u, d, s quarks we get 

q„q v P^(q) = -P L (q 2 )q 2 = N 2 (a s /4n) 2 i J d^xe^x 

x (o\tg; u (x)g; u (x),gz(o)gz(oM> (?) 

where G" is the gluonic field strength, G^ v = (l/2)e t _ llJ x a Gx a . At q 2 — ► 0, we 
have 

f 2 = -(1/3)P L (0) = ^N 2 fX '(0), (8) 
where x(q 2 ) is the topological susceptibility 

X {q 2 )=i I d A xe^{Q\TQ b {x),Qs{m) (9) 



Q 5 (x) = (a s /8n) G^(x)G^(0), (10) 

x(0) = if_thcre is at least one massless quark i. 

In ref. □ the sum rule, expressing S in terms of /g (4) or x'(0) was found. 
The OPE up to dimension d = 7 was performed. 



However, the accuracy of the calculation was not good enough for reliable 
calculation of E in terms of f§: the necessary requirement of the method - the 
weak dependence of the result on the Borel parameter was not well satisfied. 

In this paper we improve the accuracy of the calculation by going to higher 
order terms in OPE up to dimension 9 operators. Under the assumption of 
factorization - the saturation of the product of four-quark operators by the 
contribution of an intermediate vacuum state - the dimension 8 v.e.v.'s are 
accounted (times A^): 

- q(0 | qa a0 (l/2)X n G n a0 q ■ qq \ 0) = m§<0 | qq | 0) 2 (11) 

where m 2 , = 0.8 ± 0.2 GeV 2 was determined ini The sum rule for E is given 
by 

E + C Q M 2 = -1 + ^e m2 / A/2 |a 2 L 4 / 9 + +6tt 2 f 2 M i E l (j^j iT 4 / 9 + 

Here M 2 is the Borel parameter, A at is defined as A^ = 327r 4 A 2 v — 
2.1 GeV 6 , {0 1 77 1^3} = Xnv p , where v p is proton .spinor, W 2 is the continuum 
threshold, W 2 = 2.5 GeV 2 , h = 3 • 10" 4 GeV 4 E, 

a = — (27r) 2 {0|gg|0) = 0.55 GeV 3 (13) 

E (x) = l-e- x , E 1 {x) = l-{l + x)e- x 

L = ln(M/ A)/ln(fi/ A), A = Aqcd = 200 MeV and the normalization point 
fi was chosen /i = 1 GeV. When deriving (12) the sum rule for the nucleon 
mass was exploited what results in appearance of the first term, -1, in the right 
hand side (rhs) of (12). This term absorbs the contributions of the bare loop, 
gluonic condensate as well as a s corrections to them and essential part of terms, 
proportional to a 2 and m§a 2 . The values of the parameters, a,\ N , W 2 taken, 
above were chosen by the best fit of the sum rules for the nucleon mass (seeB, 
Appendix B) performed at A = 200 MeV. The unknown constant Co in the 
left-hand side of (12) corresponds to the contribution of inelastic transitions 
p — > N* — > interaction withA^ — + p (and in inverse order). The sum rule (12) 
as well as the sum rule for the nucleon mass is reliable in the interval of the 
Borel parameter M 2 0.85 < M 2 < 1.4 GeV 2 . The Af 2 -dependence of the rhs of 
(12) at f§ = 3 x 10~ 2 GeV 2 is plotted in Fig. la. The complicated expression m 



rhs of (12) is indeed an almost linear function of M 2 in the given interval! The 
best values of E = E /l * and C = C$ U are found from the \ 2 fitting procedure 



X 2 = - J2 - C o U Mf R(M 2 )] 2 = mm, (14) 

i=i 

where i?(M 2 ) is the rhs of (12). 

The values of E as a function of f§ are plotted in Fig. lb together with yx^. 
In our approach the gluonic contribution cannot be [Separated and is included 
in E. The experimental value of E can be estimated&B as E = 0.3 ± 0.1. Then 
from Fig.lb we have f 2 = (2.8 ± 0.7) x 10~ 2 GeV 2 and x'(0) = (2.3 ± 0.6) x 
10~ 3 GeV 2 . The error in /q and \' besides the experimentall error includes 
the uncertainty in the sum rule estimated as equal to the contribution of the 
last term in OPE (two last terms in Eq.12) and a possible role of NLO a s 
corrections. 

From x 2 fit, i-e. from the requirement of the selfconsistancy of the sum 
rule, allowing the deviation of x 2 by a factor 1.5 from its minimal value, we 
find: E = 0.33 ± 0.20, x'(Q) = (2.4 ± 0.5) • 10~ 3 GeV 2 . 

From the sum rule analogous (12) it is possible to find g\ - the proton 
coupling constant with the octet axial current, which enters the QCD formula 
for T p , n . 

The M 2 -dependence of g\ + C%M 2 is presented in Fig. la and the best fit 
according to the fitting procedure (14) at 1.0 < M 2 < 1.3 GeV 2 gives 

g\ = 0.65 ±0.15, C 8 = 0.10 GeV~ 2 y/x 2 = 1.2 x 10~ 3 (15) 

The obtained value of g\ within the errors coincides with g\ = 0.59 ± 0.02 
found from the data, on baryon octet /3-decays under assumption of strict SU(3) 
flavour symmetry a. 

A similar sum rule with the account of dimension 9 operators can be 
derived also for gA - the nucleon axial /3-decay coupling constant. It is an 
improvement of the sum rule found inE3 and has the form 



g A +C A M 2 = l+° e ™ 2 / M2 
9X 2 N 



2 7-4/9 , o 2 2.f2 7i /f 2 1 2 TO ^ f 2 ^ 

a L +2ir mJ^M --a jj^a + 3™'/^ ]jf2 



The M 2 dependence of g A — 1 + C A M 2 is plotted in Fig. la, lower curve. The 
best fit gives 



g A = 1.37 ± 0.10, C A = -0.088 GeV~ 2 , y/x 2 = 1.0 x 10~ 3 (17) 




Figure 1: a.) The M 2 -dependence of E + C M 2 at / 2 = 3 X 10" 2 GeV 2 , g\ + C 8 M 2 , and 
gA — 1 + C^M 2 . b.)S (solid line) and -^/x 2 (crossed line) as a functions of /q . 



in comparison with the world average gA — 1.260 ± 0.002. 
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